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The decay of 56 Ni* , formed in 32 S + 24 Mg reaction at the incident energies E cm =5l.Q and 60.5 
MeV, is calculated as a cluster decay process within the Preformed Cluster-decay Model (PCM) of 
Gupta et al. re-formulated for hot compound systems. Interesting enough, the cluster decay process 
is shown to contain the complete structure of both the measured fragment cross sections and total 
kinetic energies (TKEs). The observed deformed shapes of the exit channel fragments are simulated 
by introducing the neck-length parameter at the scission configuration, which nearly coincides the 
56 'Ni saddle configuration. This is the only parameter of the model, which though is also defined 
in terms of the binding energy of the hot compound system and the ground-state binding energies 
of the various emitted fragments. For the temperature effects included in shell corrections only, the 
normalized Q-nucleus s-wave cross sections calculated for nuclear shapes with outgoing fragments 
separated within nuclear proximity limit (here ~0.3 fm) can be compared with the experimental 
data, and the TKEs are found to be in reasonably good agreement with experiments for the angular 
momentum effects added in the sticking limit for the moment of inertia. The incident energy effects 
are also shown in predicting different separation distances and angular momentum values for the 
best fit. Also, some light particle production (other than the evaporation residue, not treated here) 
is predicted at these energies and, interestingly, 4 He, which belongs to evaporation residue, is found 
missing as a dynamical cluster-decay fragment. Similar results are obtained for temperature effects 
included in all the terms of the potential energy. The non-a fragments are now equally important 
and hence present a more realistic situation with respect to experiments. 
PACS Nos. 25.70.Jj, 23.70. +j, 24.10.-i, 23.60.+e 



I. INTRODUCTION 



Experimentally, 5e Ni is an extensively studied compound system by using different entrance channels, namely 
16q _|_40 (j a ^ 28 _|_28 anc j 32g _|_24 anc j var i ous incident energies ranging from 1.5 to 2.2 times the Coulomb 
barrier (see, e.g., the review [1] and the other direct and more recent Refs. [2-8]). At such incident energies, the 
incident flux is found to get trapped by the formation of a compound nucleus (CN), which is in addition to a 
significant large- angle elastic scattering cross section. For lighter masses (Aqn < 44), such a compound nucleus 
decays subsequently by the emission of mainly light particles (n, p, a) and 7-rays; i.e. with very small component of 
heavy fragment (A > 4) emission. An experimental measure of this so-called particle evaporation residue yield is the 
CN fusion cross section. For somewhat heavier systems, like 48 Cr and 56 Ni, a significant decay strength to A > 4 
fragments, the mass-asymmetric channels, is also observed which could apparently not arise from a direct reaction 
mechanism because of the large mass- asymmetry differences between the entrance and exit channels. The measured 
angular distributions and energy spectra are consistent with fission-like decays of the respective compound systems. 

For the 32 S + 24 Mg ^ 56 Ni* reaction, in one of the experiments, the mass spectra for A=12 to 28 fragments and 
the total kinetic energy (TKE) for only the most favoured (enhanced yields) a-nucleus fragments are measured at 
the energies Ei a b = 121.1 and 141.8 MeV, or equivalently at E cm =51.6 and 60.5 MeV, respectively [2,3]. Note that 
56 Ni is a negative Q-value system (negative Q ou t, different for different exit channels) and hence would decay only if 
it were produced in heavy ion reactions with sufficient compound nucleus excitation energy Eq N (= E cm + Qi n ), to 
compensate for the negative Q ou t, the deformation energy of fragments Ed, their total kinetic energy TKE and the 
total excitation energy TXE, in the exit channel, as 

E* CN =1 Qout{T) I +E d (T) + TKE(T) + TXE(T); (1) 

see Fig.l where Ed is neglected because fragments are considered to be spherical. Here Qi n is positive (=16.68 MeV 
for 32 S + 24 Mg entrance channel) and hence adds to the entrance channel kinetic energy E cm of the two incoming 
nuclei in their ground states. In another experiment [4] for 32 S + 24 Mg reaction at E cm =51.0 and 54.5 MeV, the 
excitation-energy spectra for only the symmetric 2S Si + 28 Si and near-symmetric Mg + 32 S channels are measured, 
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whose analysis indicate that a specific set of states in 2S Si correspond to highly deformed bands. In other words, the 
expected shapes of some of the observed fragments in the exit channel could be relatively deformed. It is interesting 
to note that this result is supported by a very recent study of the 28 Si + 28 Si reaction at E cm = 55 MeV, where the 
population of highly excited states in the 24 Mg, 28 Si and 32 S nuclei indicated a selective and enhanced population of 
deformed bands [6]. In a still other recent experiment [7], the incident energy used in the same 32 S + 24 Mg ^ 56 Ni* 
reaction is Ei a b = 130 MeV and an enhanced emission yield by a factor of 1.5 to 1.8 is observed for 8 Be over the two 
a-particles. The aim of our present work is to understand some of the results of these experiments. 

The above stated light particles (A < 4) production, the evaporation residue, is very satisfactorily understood as 
the cquliberated compound nucleus emission in the statistical Hauser Feshbach analysis [3,9-12], using the LILITA 
or CASCADE codes. The Hauser Feshbach calculations are also extended to include the complex fragments, like the 
ones observed in the experiments mentioned above. These are considered in the, so-called, BUSCO code [10] or the 
Extended Hauser- Feshbach scission-point model [12]. Within the framework of the Extended Hauser- Feshbach method 
[12], the above noted observed enhanced emission of s Be over the evaporation of two a-particles in 32 S+ 24 Mg reaction 
is shown related to an increased deformation of the heavier fragment 48 Cr [7]. The emission of complex fragments 
(A > 4, also called the intermediate mass fragments, IMFs, or "clusters") is alternatively treated as the binary 
fission of a compound nucleus in the statistical fission models [13,14], using the GEMINI code [9] or the saddle-point 
"transition-state" model [3,5,11]. The transition-state model, treating the complex fragments emission as a compound- 
nucleus fission process (the fusion-fission) seems to explain the observed mass spectra and excitation-energy spectra 
rather well for the 32 S + 24 Mg reaction at the two energies used in respective experiments [3,4]. Also, the measured 
TKE for the symmetric fission is comparable to the saddle-point potential energy at I = 36?i [3]. Then, there are 
other processes, like the deep-inelastic (DI) orbiting or scattering, that have also been studied for this reaction but 
do not seem to explain the observed data [3] . 

In the statistical fission models [13,14], the fission decay of a compound nucleus is determined by the phase space 
(level density) available at the "transition" configuration, which is saddle or scission in these models. For light 
systems, this choice can lead to a significant population of many energetically allowed mass channels, though there 
is no structure information of the compound system in these fission models. However, the structure effects of the 
compound system seem to influence the observed yields strongly since strong resonance behaviour is observed in the 
measured excitation functions of large-angle elastic and inelastic scattering yields in several light systems (see, e.g., 
[6] ) . One possibility to account for such structure effects is via the process of fragments (or clusters) preformation in a 
compound nucleus and its subsequent decay as a cluster decay process, proposed recently by some of us [15,16]. The 
structure information enters the process via the preformation probability (also, known as the spectroscopic factors) 
of the fragments. We follow this approach of preformed cluster decay [15,16] here in this paper. 

The cluster decay process was recently studied [15] for the compound system 56 Ni*, using the preformed cluster- 
decay model (PCM) of Gupta and collaborators [17-21]. It was shown that for the decay of 56 Ni* , the two processes 
of binary fission (the dynamical collective mass transfer calculated, by some of us [22-25], in the quantum mechanical 
fragmentation theory [26-28]) and cluster decay are almost indistinguishable, particularly at higher angular momenta. 
However, this work was a simple model study where the role of TKE was analysed and found to be significant for 
a-nucleus structure in the measured yields. This model is more recently re-formulated [16] for the IMFs emitted from 
an excited 116 Ba* compound nucleus produced in low energy 58 Ni + 58 Ni reaction. The IMFs in lle Ba* are shown 
to be produced as multiple "clusters" of masses A < 20 and only at Ei a b > 200 MeV, in agreement with experiments. 
Both of these works [15,16] show that the IMFs in decay of excited 116 Ba* or the complete mass spectra in decay of 
excited 5e Ni* have their origin in the macroscopic liquid drop energy (the shell effects are almost zero at the excitation 
energies involved). For 116 Ba* decay, the light particles (Z < 2) emission, other than the promptly emitted via the 
statistical evaporation process (not treated in this model), is also shown to be given, but at higher energies where 
only the pure liquid drop energies enter the calculations. Thus, the macroscopic liquid drop energy (Vldm) is shown 
playing the most important role in the cluster decay calculations. Apparently, the compound nucleus being hot at 
the energies involved, the Vldm should also depend on the temperature T. This is done here in this paper for the 
decay of 56 Ni* formed in 32 S + 24 Mg reaction at the two energies, E cm =51.6 and 60.5 MeV [2,3]. Also, the other 
terms of the potential, that constitute the scattering potential V(R), are considered T-dependent. 

The T-depcndcnt liquid drop model used is that of Davidson et al. [29] which is based on the semi-empirical mass 
formula of Seeger [30] . The model parameters of Seeger's formula at T=0 are re- fitted in view of the present availability 
of a larger data set for binding energies [31]. For the T-dependence in V(R), we follow Davidson et al. [29] and some 
other authors [32], discussed below. The deformation effects of the fragments (and the neck formation between them) 
are included here within the extended model of Gupta and collaborators [33-35] , via a neck- length parameter at the 
scission configuration which simulates the two centre nuclear shape parametrization, used for both the light and heavy 
nuclear systems. A similar method has been used earlier by other authors [3,11,12], discussed below. 

The dynamical cluster decay model for hot compound systems, a re-formulation of the preformed cluster-decay 
model (PCM) of Gupta and co-workers [17-21] for ground-state decays, is presented in section 2 and its application 
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to the hot 56 Ni* nucleus data from Refs. [2,3] in section 3. The (statistical) evaporation of light particles, that 
occur promptly before the beginning of the binary decay process of cluster emission studied here, is not included in 
this paper. Hence, any discussion of light particles emission is that of one which is in addition to the ones emitted 
promptly. Finally, a summary of our results is presented in section 4. 



II. THE DYNAMICAL CLUSTER DECAY MODEL FOR HOT COMPOUND SYSTEMS 

The cluster decay model developed here is the preformed cluster-decay model (PCM) of Gupta ct al. [17-21] for 
the ground-state decays, re-formulated for hot and excited compound systems. In this model, we treat the complex 
fragments (the IMFs or clusters) as dynamical collective mass motion of preformed fragments through the barrier. It 
is based on the well known dynamical (or quantum mechanical) fragmentation theory [26-28] developed for fission and 
heavy ion reactions, and used later for predicting the exotic cluster radioactivity [36-38] also. This theory is worked 
out in terms of the collective coordinates of mass asymmetry r\ = (Ai — A 2 )/(Ai + A 2 ) and relative separation R, 
which in a PCM allows to define the decay half-life Ti , or the decay constant A, as 

A = ^ = W, (2) 

2 

where Po, the preformation probability, refers to 77- motion and P, the penetrability, to R-motion. Apparentlty, the 
two motions are taken as decoupled, an assumption justified in our earlier works [26,27,39]. The f is the barrier 
assault frequency. In terms of the partial waves, the decay cross section 

a=^p^+l)P Q P- * = ^f^ (3) 

with ll = [AiA 2 /{Ai + A 2 )]m = jAm(l — ri 2 ) as the reduced mass and £ c , the critical (maximum) angular momentum, 
defined later, m is the nucleon mass. This means that A in (2) gives the s-wave cross section, with a normalization 
constant instead of the 7r//c 2 in (3). However, in the present calculations, made for £ — case, the normalization 
constant is obtained empirically from the experimental data. 

For 77-motion, we solve the stationary Schrodinger equation in 77, at a fixed R, 

'"i^l;^^*™™-*"*™ (4 » 

with v=0,l,2,3... and R = R a = Ct(= C\ + C 2 ), the first turning point, fixed empirically for the ground-state (T=0) 
decay since this value of R (instead of the compound nucleus radius Rq) assimilates to a good extent the effects of both 
the deformations Pi of two fragments and neck formation between them [35]. In other words, the deformation effects 
of the two fragments are included here in the scattering potential V(R,T=0) for each r\ by raising the first turning 
point R a from R a = Rq to R a = Ct or Ct + X) SR(/3i), which is equivalent of lowering of the barrier, as is found to be 
the case for deformed fragments [35] . This method of inclusion of fragment deformation and the paramctrization of 
the neck zone via a neck-length parameter 5R in the present calculations is quite similar to what has been achieved in 
both the transition-state model of Sanders [3,11] (in saddle point configuration) and the Extended Hauser-Feshbach 
Method of Matsuse and collaborators [12] (in scission point configuration). It is also shown in [35] that the alternative 
of calculating the fragmentation potential V(rj) and scattering potential V(R) for deformed nuclei is not practical 
since the experimental deformation parameters for all the possible fragments (Ai,A 2 ), required for calculating V(r]), 
are generally not available. The deformation effects of nuclei in our calculations are further included via the Sussmann 
central radii Ci — Ri — (b/Ri), with the radii Ri = 1.28Ay 3 — 0.76 + 0.8^ 1 ^ 3 fm and surface thickness parameter 
b =0.99 fm. Note that the C t are different for different i]- values and hence C t is C t {rj). 
The eigen-solutions of Eq. (4) give the preformation probability 

|2 



B m I I (2M) , (5) 

(i=l or 2), where tfj(rj) is V' ly=0 ('7) if the ground-state solution is chosen. However, the decay of 56 Ni in the ground-state 
(T=0, R a = C t ) is not allowed since Q ou t(T = 0) is negative. 

For the decay of a hot compound nucleus, we use an ansatz [16] for the first turning point, 

R a = C t ( V ,T) + AR( v ,T), (6) 
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which depends on the total kinetic energy TKE(T). The corresponding potential V(R a ) acts like an effective, positive 
Q-value, Qeff, for the decay of the hot compound system at temperature T to two fragments in the exit channel 
observed in the ground-states (T=0). Thus, in terms of the respective binding energies B, Q e ff is defined as 

Q eff (T) = B(T) - [B^T = 0) + B 2 (T = 0)] 
= TKE(T) 

= V(R a ). (7) 

Since, R a = Ct{rf) for T=0, AR{rf) corresponds to the change in TKE at T with respect to its value at T — and 
hence can be estimated exactly for the temperature effects included in the scattering potential V(R). Note that in Eq. 
(6) C t is also taken to depend on temperature, as is defined later in the following. Also, AR depends on rj. In the 
following, however, based on our earlier work [16], instead, we use a constant average value AR, independent of rj, 
which also takes care of the additional effects of the deformations of fragments and neck formation between 

them. Note that AR is the only parameter of the model, though it is shown that the structure of the calculated mass 
spectrum is nearly independent of the exact choice of this parameter value. The corresponding Q e ff is denoted as 

Qeff(AR). 

In the above definition of Q e ff, apparently the two fragments would come out of the barrier and go to ground state 
(T — ► 0) only by emitting some light particle(s) and/or 7-rays of energy, defined as (see Fig. 8) 



Eq. (8) means that one can also write 



E x =B{T)- B(0) 

= Qout(T)-Q out (T = 0) + AB 
= Qeff(T)-Q ma (T = 0) 
= TKE(T) - TKE(T = 0). 



(8) 



Q eff (T) = TKE(T) 

= Qout(T = 0) + E x 

= TKE(T = 0)+E x , (9) 

which is what one observes experimentally i.e. the fragments in the ground state with Q ou t{T = 0) {=TKE{T = 0)) 
and light particle(s) and 7-rays of energy E x . The remaining excitation energy of the decaying system is then, 

E* CN - E x =| Q out (T) I +TKE(T = 0) + TXE(T), (10) 

which again shows that the exit channel fragments are obtained with their TKE in the ground-state, i.e. with 
TKE(T=0). The excitation energy TXE(T) in (10) is used in the secondary emission of light particles from the frag- 
ments, which are not treated here. Instead, we compare our calculations with the primary pre-secondary-evaporation 
fragments emission data. 

We notice from Eq. (7) that for the ground-state (T=0) decay, 

Qeff(T = 0) = Q out (T = 0) = TKE(T = 0), (11) 

as is the case for exotic cluster radioactivity [21,38]. In fact, one can write Eq. (7) as 

Q eff (T)=Q out (T) + AB, (12) 

where 

AB = [Bi(T) + B 2 (T)] - [Bi(T = 0) + B 2 (T = 0)], (13) 

the difference of binding energies at temperature T and the ground-state binding energies of the two fragments. Also, 
for the ground-state (T=0) decays, according to Eq. (8), E x = (no particle or 7-ray emission), as is known to be 
true for exotic cluster radioactivity [21,38]. 

Thus, at temperature T, the preformation factor P in Eq. (5) is calculated at R a — C t (rj) + AR, with the 
temperature effects also included in tp(r]) through a Boltzmann-like function 

00 

I V> | 2 = E I ^ I' expi-E^/T), (14) 

u=0 
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with the compound nucleus temperature T (in MeV) related as 

E* CN = (A/9)T 2 -T; (15) 

and for the penetrability P, Eqs. (6) and (7) for each r\ and T-values, mean that 

V(R a ) = V{C t + AR) = V(R h ) = Q ef f(AR) = TKE(T), (16) 

with Rb as the second turning point, and penetrability P calculated as the WKB tunnelling probability for the path 
shown in Fig. 1 (or Fig. 8), as 

P = exp[~ {MV{R) - Qe f f}} 1/2 dR}, (17) 

solved analytically [18]. 

The fragmentation potential Vp(r],T) at any temperature T, in Eq.(4), is calculated within the Strutinsky renor- 
malization procedure, as 

V R ( V ,T) = J2[VLDM(A i ,Z i ,T)] +J2[sUi]ex P (-^) 
i=i i=i 

+ E c (T) + V P (T) + Ve(T), (18) 

where the T-dependent liquid drop energy Vldm(T) is that of Ref. [29] with the (Seeger's) constants at T=0 re-fitted 
to give the experimental binding energies B [31], defined as B — Vldm(T = 0) + SU. The shell corrections 5U are 
calculated in the "empirical method" of Myers and Swiatecki [40]. Some of these details are given in Appendix I. 
Figure 2 illustrates the kind of comparisons obtained for V(rj) calculated at R = C\ + C2 = Ct and T=0 for the 
experimental and newly fitted binding energies. Apparently, the binding energies fit within 1 to 1.5 MeV. 

The Vp is an additional attraction due to the nuclear proximity potential [41], which is also considered temperature- 
dependent here, 

V P (R, T) = 47ri?(T) 7 6(T)$( s , T), (19) 

where R(T) and $(s,T) are, respectively, the inverse of the root mean square radius of the Gaussian curvature and 
the universal function which is independent of the geometery of the system, given by 



a , T , f -i(s- 2.54) 2 - 0.0852(s-2.54) 3 for s < 1.2511 
[S ' ' ~ \ -3.437ea;p(-( J %) for s > 1.2511 



(20) 



R(T) 

and 7 is the specific nuclear surface tension given by 



d(T)C 2 (T) 
C t (T) ' 



7 = 0.9517 



1 - 1.7826 



N — Z 



MeVfrn 



(21) 



(22) 



In Eq. (20), s(T) (= — &(rj T ^ ) ^ s over l a P distance, in units of b, between the colliding surfaces. The temperature 
dependence in radii Ri is given as [29,32], 

R t {T) = r (T)Af = 1.07(1 + 0.01T)AJ (23) 

with the surface width 



b(T) = 0.99(1 + 0.009T 2 ). 



(24) 



The same temperature dependence of R(T) is also used for Coulomb potential E C {T) — ZiZ 2 e 2 /i?(T), where 
the charges Zi are fixed by minimizing the potential Vp(i],T) in the charge asymmetry coordinate r\z = 



5 



[Z\ — Z 2 )/(Zi + Z 2 ). The shell corrections SU in Eq. (18) are considered to vanish exponentially for T = 1.5 
MeV [42]. 

Also, for the angular momentum effects (so far included here for the calculation of total kinetic energy only) 

In the non-sticking limit, where R a — C\{T) + C 2 (T) + AR = C t (T) + AR, the moment of inertia in (25) is given by 

I(T) = I NS (T) = nR 2 a . (26) 

In this case, the separation distance AR is assumed to be beyond the range of nuclear proximity forces, which is about 
2 fm . However, when AR is within the range of nuclear proximity (< 2 fm), we get in the complete sticking limit 

I(T) = I S (T) = vR 2 a + 2 -A x mCl + 2 -A 2 mC\. (27) 
For the I- value, in terms of the bombarding energy E cm of the entrance channel we have 

e = t c = Ra VMEcm-V(R a , mn J. = 0)]/h, (28) 

or, alternatively, it could be fixed for the vanishing of fusion barrier. In this work, however, we use I = for the IMF 
cross sections and take £ c as a variable parameter for total kinetic energy (TKE) calculations (see Fig. 8). 

The mass parameters B vv (ri), representing the kinetic energy part in Eq. (4), are the smooth classical hydrody- 
namical masses [43], since we are dealing here with a situation where the shell effects are almost completely washed 
out. 

The assault frequency v , in Eq.(2), is given simply as 

1/0 R ' (29) 

with the kinetic energy of the lighter fragment E 2 — (Ai/A)Q e ff, for the Q e ff shared between the two fragments as 
inverse of their masses. However, for the calculations of s-wave cross sections, instead of we use an empirically 
determined normalization constant. 

Finally, the temperature dependent scattering potential V(R,T), normalised to the exit channel binding energy, is 

V(R,T) = Z!Z 2 e 2 /R{T) + V P (T) + V e (T). (30) 

This means that all energies are measured w.r.t B\(T) + B 2 (T), and the fragments go to ground state (T — ► 0) via 
the emission of light particle (s) and/or 7-rays of energy E x . 



III. CALCULATIONS 



The calculations are made in two steps: (i) with temperature effects included only in the shell corrections, i.e., 
using SU(T) 7 but T-independent Vldm and V(R); and (ii) with temperature effects included also in both the liquid 
drop energy and scattering potential, i.e., using Vldm(T), 8U(T) and V(R,T). This allows us to study explicitly the 
role of temperature in different terms of the potential. In both sets of the calculations, we first take £=0, i.e., use 
V(=0 through out, but then study the effect of adding this term to the potential V(R) for calculating the total kinetic 
energy (TKE) alone. 



A. Temperature effects only in shell corrections 

Figure 3 gives our calculated fragmentation potentials V{q, T) for 56 Ni* at T=0, as well as at other two temperatures 
refering to the compound nucleus excitation energies Eq N of the experiments of Ref. [2,3]. The R- values chosen are: 
R = R a = C t a,t T=0, and, as before [16], R = R a +~AR with ~AR =0.30 and 0.31 fm, arbitrarily, for T =3.39 and 3.60 
MeV which correspond to the experimental energies E cm =51.6 and 60.5 MeV, respectively. The near independence 
of the structure in V(rj) on R-value was studied in our earlier works [21,28]. The SU at these temperatures reduce 
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almost to zero. However, we notice that the N=Z, A=4n a-nucleus structure is obtained at all temperatures, which 
has its origin apparently in the macroscopic liquid drop energy and is due to the "Wigner term" in it, as was also 
shown earlier in Refs. [15,16]. Note that here the Vldm and other R-depcndcnt terms (E c and Vp) are not yet 
T-dependent (see next subsection). This means that for use of only SU as T-dependcnt, the N=Z a-nuclci fragments 
should be produced preferentially in the decay of 56 Ni* at all temperatures. 

The preformation probability Po of the fragments, calculated for the potentials in Fig. 3, is given in Fig. 4. The 
case of T = is not shown here since cold 56 Ni (in the ground-state) can not decay because of its negative Q-value. 
Interesting enough, for both the temperatures (the two temperatures are nearly the same), the yields are large for 
only a small window of A < 16 fragments, including the light particles (A < 3). Also, the a-nucleus fragments 
4 He, 8 Be, 12 C and 16 and the light particle 1 H (in addition to the evaporation residues, not included here) are 
preferentially preformed. This means that, out of all the fragments observed in the decay of 56 Ni*, the ones with 
A < 16 arc strongly preformed. The other ones with A > 16, if observed, must have larger penetrability P, since the 
decay constant is a combined effect of both the preformation factor Po & n d penetrability P (vq is nearly constant). 

Figure 5 gives the results of our calculation for the normalized decay constants, equivalently, the s-wave production 
cross sections for only the most favoured (largest yields or cross sections) a-nucleus fragments, compared with the 
experimental data at two energies, taken from Fig. 9 of Ref. [3]. In the lower panel, the calculation at E cm —51.6 MeV, 
using AP=0.3 fm, is fully normalized to the experimental data for the favoured a-nucleus fragments only. Then, in 
the upper panel, for the higher energy P cm =60.5 MeV we find that, for the use of the same normalization as obtained 
in lower panel and for a further normalization of the A=12 fragment yield, the best fit to the a-nuclei fragment data is 
obtained for AP=0.29 fm, a value lower than that used for the lower incident energy P cm =51.6 MeV. This is contrary 
to the expected behaviour of increased R at higher temeratures but, as we shall see below in Fig. 7, this is a result of 
our having not included here the contribution of angular momentum term in the fragmentation potential (Vg = in 
V(r],T)) and hence in the cross sections. Also, the inclusion of temperature effects in other terms (the Vldm, E c and 
Vp) are important, as is shown below in subsection III.B. Hence, Fig. 5 (and Fig. 7 below) show that the dynamical 
cluster-decay model contains the required structure of the measured yields (and TKEs) in this experiment [2,3]. 

Figure 6 shows the complete mass spectra for decay of 56 Ni* calculated at both the energies and compared with 
the measured yields [3]. The calculated yields are for the energetically favoured, most probable, mass fragments (see 
Figs. 3 and 4). Note that the experimental data in [2,3] are available only for fragments heavier than mass 11, and in 
steps of mass one for P cm =60.5 MeV, but in steps of only mass two for P cm =51.6 MeV due to a deteriorated mass 
resolution at the lower bombarding energy. For comparisons, the calculations are normalized to the experimental data 
for one fragment mass (A=20) only. The role of the penetrability P is evident in this figure, since some of the strongly 
preformed fragments, like 4 He and 1 H in Fig. 4, are now shown as less favoured decays (smaller cross sections, not 
shown in Fig. 6 since they lie below the chosen scale). The same is true of weakly preformed fragments (in Fig. 4), 
with A > 16. Specifically, amongst the light particles, mass 3 fragment ( 3 He) is shown to be produced with a large 
cross section, and for lighter fragments (A < 12), instead of A=8 ( 8 Pe), the fragments with A=6 and 10 are shown 
to be produced with larger cross sections. This means that of all the residue products (A < 4, not studied here) 
only mass 3 fragment ( 3 He) is produced and that the mass 4 ( 4 He) fragment is not at all produced as a dynamical 
cluster-decay fragment. This non-occurence of 4 He as a dynamical cluster-decay product, is an interesting result, 
giving a strong support to the credential of the model. For mass 8 (Be) decay, perhaps the contribution of higher 
I- values is important. For the heavier fragments (A > 20), the calculated cross sections are rather small due to the 
fact that here the contribution of only I = term is considered. Also, in experiments it is difficult to separate the 
contributions of direct (such as alpha-transfer and orbiting processes) and compound nucleus yields for the heavy 
mass fragments (A> 20) (see Ref. [1] and references therein). Thus, in view of the fact that we are dealing here with 
only the I = case and that the temperature effects are not included in full in the potential, the comparisons in Fig. 
6 between the theory and experiments could be said atleast reasonable. 

Figure 7 shows the results of our calculation for total kinetic energy (TKE), with angular momentum £ effects 
included only in the scattering potential V(R). We notice that the calculated TKEs for the sticking limit (using Is) 
compare reasonably good with the experimental data. This means that, even though AR is non-zero (=0.29 and 0.3 
fm), the sticking limit for the moment of inertia is preferred. Also, unlike the AR- values, the I- values required for the 
case of higher energy data is now of a larger value (£ = 25h for E cm =60.5 MeV as compared to 24?i for E cm =51.6 
MeV), as expected. The measured TKEs are taken from Ref. [3]. 

B. Fully temperature-dependent potential 

Figure 8 shows the scattering potential V(R,T,£) for temperature effects included in all the terms of the potential 
(compare this figure for I — 0, with Fig. 1 where temperature effects are included in SU only). Notice that as lvalue 
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increases, the TKE(Ai?)-value increases, since the decay path for all the t- values begins at R = R a . Figure 9 gives 
our calculated fragmentation potentials V{rj^T). The T- values chosen are the same as in Fig. 3, where temperature 
effects were included only in the shell corrections. The R- values here are: R{T) — Ci(T) + C2(T) + AR = C't(T) + AR, 
with Ai?-values as shown in the figure. We notice in Fig. 9 that, due to the inclusion of temperature effects in all 
terms, the minima in the potential, which were earlier only for a-nuclei, are now obtained for both the a and non-a 
fragments. This happens, possibly, due to the pairing energy term S(T) in formula (32) of Davidson et al. [29], which 
goes to zero for T > 2 MeV. Thus, with the addition of temperature, not only the shell structure effects go to zero but 
also the explicitly preferred a-nucleus structure washes out. Also, we notice that the light particles (A < 4) structure 
changes; in particular, the minimum at A He disappears and a shallow minimum at 2 H appears. 

Figure 10 gives the preformation factors Pq for the two experimentally chosen temperatures only, since the ground- 
state (T=0) decay is not possible. We notice that the formation yields are large only for light fragments (A < 16) 
and are of the same orders as in Fig. 4, except that now the non-a fragments are also preformed equally strongly. 
However, the calculated decay constants, equivalently, the fragments (s-wave) production cross sections, in Fig. 11 do 
not show much improvement in their comparisons with experiments. The comparisons are now somewhat better for 
the heavier fragments but the yields for fragments lighter than A=9 are very low, lying below the chosen scale. On 
the other hand, the calculated TKEs in Fig. 12 compare nicely (even better than in Fig. 7) with the experimental 
data. Only the case of sticking limit is shown since the Ai?-values are still within the proximity limits. Note that the 
^-dependent contribution is so far added here only in the scattering potential V(R, T) and not yet in the fragmentation 
potential V(r], T) which is needed for both the preformation factor and penetrability. This extention is being carried 
out. 



IV. SUMMARY 

In summary, we have reformulated for hot nuclear systems, the preformed cluster-decay model (PCM) of Gupta 
and collaborators for ground-state decays and applied it for the first time to the decay of a light compound nucleus 
such as 56 Ni* formed in the reaction 32 S + Mg carried out at two incident energies E cm =51.6 and 60.5 MeV [2,3]. 
In this experiment, the mass spectra for fragments heavier than mass 12 and the total kinetic energies (TKEs) for 
only the favoured a-nucleus fragments are measured. Also, at another energy, in between the two above, an enhanced 
yield is observed for 8 Be over the two a-particles emission [7]. Our calculations are made first for the temperature 
effects included only in shell corrections and then in all terms of the potential, and in each case for 1=0 only. The 
contribution due to I is added only for estimating the TKEs. Similar to the saddle-point model [11] and/ or the 
scission-point model [12], the deformations of the fragments are taken into account by the parametrization of the 
neck- in zone, proposed by Gupta and collaborators [33-35] . This quantity is ^-dependent and could be calculated but 
is taken as a parameter here, which is the only parameter of the model. 

For the temperature effects included in shell corrections only, we find that the a-nucleus fragments are favourably 
preformed and are due to the macroscopic liquid drop energy alone since the shell effects are almost zero at the 
energies under consideration. The calculated decay constants or the normalized s-wave cross sections, in particular 
for the a-nucleus fragments are found to contain the complete structure of the experiments for a nuclear shape with 
fragments separated by about 0.3 fm which is within the limits of nuclear proximity effects. Some of the light particles 
(other than the ones constituting the evaporation residue, not included here) are also predicted to be there in the 
mass spectra, but 4 He is shown to be absent. With angular moments effects included, the calculated TKEs are found 
to compare rather nicely with experimental data for the moment of inertia calculated for a sticking limit. 

For the full temperature effects in the potential, the non-a fragments are also preformed equally strongly as the 
a-nucleus fragments. The cluster decay process now occurs at a somewhat larger separation distance, which is also 
temperature dependent. Hence, the TKEs for a sticking moment of inertia are now in somewhat better agreement 
with the experiments. However, the comparison between the calculated (s-wave) and measured mass spectra is not 
improved much, which calls for the inclusion of ^-dependent potential in the calculations of yields also, which is 
underway 

V. APPENDIX I: TEMPERATURE DEPENDENT BINDING ENERGIES 

In Eq. (18) we have defined, within the Strutinsky renormalization procedure, the binding energy B of a nucleus 
at temperature T as the sum of liquid drop energy Vldm(T) and shell correction 5U(T), 

B{T) = V LDM (T) + SUexp(-^). (31) 

Jo 
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The T dependent liquid drop part of the binding energy Vldm(T) used here is that of Davidson et al. [29], based on 
the semi-empirical mass formula of Sccgcr [30] , as 



V LDM {T) = a(T)A + 0(T)A* + ( 7 (T) - 



Ai ) 

// 2 + 2|/|\ Z 2 / 0.7636 

x I a J + ^jH I 1 ~ ~zT~ 

2,9 ^ +5{T) ^ (32) 



lr (T)Ai} 2 J A 
where 

I = a a (Z-N), o =l, 
and, respectively, for even-even, even-odd, and odd-odd nuclei, 

f(Z,A) = (-1,0,1). 

For T =0, Seeger [30] obtained the constants, by fitting all even-even nuclei and 488 odd-A nuclei available at that 
time, as 

a(0) = -16.11MeV, /3(0) = 20.21MeV, 
7(0) = 20.65MeV, 77(0) = 48.00MeV, 

5(0) = 33.0MeV, 



with the pairing energy term 



from Ref. [44]. Evidently, these constants need be re-fitted since a large amount of data has become available [31], 
particularly for neutron-rich nuclei. We found that the measured binding energies could be fitted within 1 to 1.5 
MeV by changing the bulk constant a(0) and introducing a proton, neutron asymmetry constant a a . The a(0) works 
as an overall scaling factor and a a controls the curvature of the experimental parabola (and hence helps to fit the 
binding energies for neutron-rich nuclei), as expected. Table 1 gives the new a(0) and a a constants for all the known 
nuclei with 1 < Z < 28, relevant to present problem. The kind of comparisons obtained between the experimental 
and calculated binding energies is already illustrated in Fig. 2. 

The T-dependent constants in Eq. (32) were obtained numerically by Davidson et al. [29] for the available exper- 
imental information on excited states of 313 nuclei in the mass region 22 < A < 250 by determining the partition 
function Z(T) of each nucleus in the canonical ensemble and making a least squares fit of the excitation energy 

E ex (T) = V LDM (T) - V LDM (T = 0) 

to the ensemble average 

E ex (T) = T 2 -^lnZ(T). 

The a(T), (3{T), 7 (T), r](T) and S(T) thus obtained are given in Figure 1 of Ref. [29] for T < 4MeV, extrapolated 
linearly for higher temperatures. For the bulk constant a(T), instead, an empirically fitted expression to a Fermi gas 
model is used, as 

a(T) - a(0) + ^. 

Also, the S(T) is constrained to be positive definite at all temperatures, with S(T > 2MeV) =0. Finally, the analytical 
form for r (T) , taken from Ref. [45] , is 

r (T) = 1.07(1 + 0.01T). 

For the shell corrections 5U in Eq. (31), since there is no microscopic shell model known that gives the shell 
corrections for light nuclei, we use the empirical formula of Myers and Swiatecki [40]. For spherical shapes, 
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where 



su = c 



F{N) + F{Z) 
' (A/2)i 



- cA* 



(33) 



F(X) = - 



M; 



M} 



M, - Mi. 



(X-M^-^xi-Ml,) 



(34) 



with X=N or Z, M;_i < X < Mi and M { as the magic numbers 2,8,14 (or 20), 28,50,82,126 and 184 for both neutrons 
and protons. The constants C=5.8 MeV and c=0.26. In this paper, we refer to the use of magic numbers 14 or 20 as 
MS14 or MS20 parametrization. 
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FIGURE CAPTIONS 



Fig.l The s-wave {£ = 0) scattering potential for 56 Ni* -^ V2 C+ 44: Ti, calculated for no temperature effects in E c and 
Vp, i.e. V{R) = E c + Vp. The Q- values are calculated from T-dependent binding energy B(T) = Vldm + 5U{T). 
The actually calculated decay path for V(R a ) = Q e ff{AR) = V(C t (T) + AR) is shown, where AR is an average 
of the separation distances for different fragmentations (different rj- values). 

Fig. 2 The fragmentation potential for 56 Ni at T=0, R = Ct, using the experimental binding energies (solid squares) 
[31] and the empirically fitted Seeger's binding energies (solid circles) with the new constants of Table 1. Here, 
MS14 means the shell corrections from the empirical method of Myers and Swiatecki [40] with Z and N=14 as 
the magic numbers. 

Fig. 3 The fragmentation potentials V(r], R, T) for 56 Ni* compound system, calculated at the ground-state (T=0, 
R a = C t ) and at various temperatures with R a = C\{T) + C*2(T) + AR values as shown. The T-dependence is 
included only in the shell corrections. 

Fig. 4 The fragment preformation probability P for 56 Ni*, calculated by using the fragmentation potentials in Fig. 
3 for the two experimental T-values only. 

Fig. 5 The calculated s-wave cross sections for the a-nucleus fragments compared with the measured ones produced 
in the reaction 32 S + 24 Mg -^ 56 Ni* at E cm =51.6 and 60.5 MeV. The data are from Fig. 9 of Ref. [3]. The 
calculations for E cm =51.6 MeV in the lower panel are made for AR =0.30 fm and are normalized completely to 
the experimental data. Using the same normalization, the calculations for E cm =60.5 MeV in the upper panel 
are made for AR =0.29, 0.30 and 0.31 fm and compared with the experimental data, for a further normalization 
of the data at fragment mass A=12. Only the a-nucleus fragments are studied, since they have the largest cross 
sections. The dotted lines are drawn only for the guide of eyes. 

Fig. 6 Same as for Fig. 5, but studied for all the fragments at E cm =51.6 MeV, AR =0.30 fm (upper panel) and 
E cm =60.5 MeV and AR =0.29 fm (lower panel). The calculations are normalized to the experimental data for 
one fragment mass (A=20) only. The calculated (s-wave) cross sections are for the energetically most favoured 
fragments in r\ coordinate i.e. fragments lying at the minimum in the fragmentation potential V{rf), minimized 
in rjz coordinate. 

Fig. 7 The measured and calculated total kinetic energy (TKE) for average AR for the the reaction 32 5 + 24 Mg — > 56 
Ni* — > A\ + A2, at the two incident energies. The calculations for I ^ are made for both the cases of sticking 
and non-sticking limits (see text) . The data are from Fig. 5 (summed over all the angles) of Ref. [3] . The same 
data are also given in Fig. 10 of Ref. [11], where it should be noted that Fig. 10(a) refers to _E c . m .=60.5 MeV 
and Fig. 10(b) to E c , m =51.6 MeV. 

Fig. 8 Same as for Fig. 1, but with £, and T dependences included in E c and Vp also, i.e., the scattering potential 
is V{R,T,£) = E C (T) + V P {T) + V t {T) with Q-valuc now calculated from B(T) = V LDM {T) + SU(T). Only 
the sticking limit of moment of inertia is used in Ve(T). The T=0 potential is shown for comparisons. For all 
£- values, the decay path (dotted line), shown for AR, begins at R = R a (marked explicitly). The distribution 
of energies and definitions of other quantities like AB and E x are indicated for the calculated AR- value. 

Fig. 9 Same as for Fig. 3, but for T-dependence in all the terms of the fragmentation potential, and at AR values as 
shown. 

Fig. 10 Same as for Fig. 4, but for the fragmentation potential of Fig. 9. 

Fig. 11 Same as for Fig. 6, but for T-dependence in all the terms of the fragmentation potential, and at AR values 
as shown. For lighter fragments, the calculated yields are not shown as they lie below the chosen scale. 

Fig. 12 Same as for Fig. 7, but for T-dependence in all the terms of the fragmentation potential, and at AR values 
as shown. 
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Table 1. Re-fitted bulk and asymmetry constants for Seeger's mass formula. 



z 


N 


a(0) 


3-a 


z 


N 


a(0) 




Z 


N 


a(0) 




1 


2 


-15.85 


0.10 


6 


9 


-15.70 


0.10 


10 


7 


-15.70 


0.50 




3 


-16.95 


0.12 




10 


-15.10 


0.10 




8 


-15.90 


0.90 




4 


-13.00 


0.05 




11 


-14.80 


0.10 




13 


-15.95 


0.50 




5 


-13.70 


0.12 




12,13,15,16 


-15.00 


0.80 
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0.50 


2 
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0.10 




14 


-14.85 


0.80 




9-12,15-22 
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0.88 




2 


-16.00 


0.10 


7 


3 


-14.30 


0.20 


11 


7 


-15.55 


0.50 




3 


-16.80 


0.30 




4 


-15.20 


0.50 




8 


-15.80 


0.50 




4,5 


-14.20 


0.30 




5 


-16.20 


0.80 




14 


-15.95 


0.50 




6 


-13.50 


0.10 




6 


-16.55 


0.80 




9-13,15-24 


-16.20 


0.86 




7,8 
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0.10 
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Fig.1 "Cluster decay of hot 56Ni\..." - Gupta et al. 
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Fig.2 "Cluster decay of hot 56Nj\..." - Gupta et al. 
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Fig. 3 "Cluster decay of hot 56Ni\..." - Gupta et al. 
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Fig. 4 "Cluster decay of hot 56Ni\..." - Gupta et al. 
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Fig. 5 "Cluster decay of hot 56Ni\..." - Gupta et al. 
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Fig. 6 "Cluster decay of hot 56Ni\..." - Gupta et al. 
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Fig. 7 "Cluster decay of hot 56Ni\..." - Gupta et al. 
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Fig. 8 "Cluster decay of hot 56Ni\..." - Gupta et al. 



30 



> 
CD 



20- 



10- 



0- 



10- 



32S+24Mg— -^56Ni— + A +A 

1 2 



> -20- 



-30- 



-40- 



-50- 



■60- 



R =C +AR 

a t 



\ 



A-A 



a -\aa 

• 26AI 28 Si 



,-/V v^s* 



16Q 18 F 



y 



; /12 c, • 



\/ • • "14M 

• iob 



1H 2 H 



6u 8 Be 



12Q 



\ 



16Q 



-"V 

20Ne 



\ 



24Mg 



28Si 



/ 



1H 



8Be 

Refitted Seeger's B's + MS14 



4He 



— ■— T=0.000, AR=0.0 
— •— T=3.394, AR=0.41 
— T=3.604, AR=0.89 
(T in MeV; R in fm) 



T 



T 



T 







5 10 15 20 25 

Fragment mass number 



30 



Fig. 9 "Cluster decay of hot 56Ni\..." - Gupta et al. 



Fig. 10 "Cluster decay of hot 56Ni\..." - Gupta et al. 
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Fig. 1 1 "Cluster decay of hot 56Ni*..." - Gupta et al. 
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Fig.12 "Cluster decay of hot 56Ni*...." - Gupta et al. 



